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A measurement of the top-antitop production charge asymmetry Ac is presented using data corre- 
sponding to an integrated luminosity of 1 .04 fb"^ of pp collisions at v^ = 7 TeV collected by the ATLAS 
detector at the LHC. Events are selected with a single lepton (electron or muon), missing transverse 
momentum and at least four jets of which at least one jet is identified as coming from a 6-quark. A 
kinematic fit is used to reconstruct the tt event topology. After background subtraction, a Bayesian 
unfolding procedure is performed to correct for acceptance and detector effects. The measured value 
of Ac is Ac = -0.019 ±0.028 (stat.)± 0.024 (syst.), consistent with the prediction from the MC@NLO 
Monte Carlo generator of Ac = 0.006 ± 0.002. Measurements of Ac in two ranges of invariant mass 
of the top-antitop pair are also shown. 
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Abstract A measurement of the top-antitop produc- 
tion charge asymmetry Ac is presented using data cor- 
responding to an integrated luminosity of 1.04 fb~^ of 
pp collisions at y/s = 7 TeV collected by the ATLAS 
detector at the LHC. Events are selected with a single 
lepton (electron or muon), missing transverse momen- 
tum and at least four jets of which at least one jet is 
identified as coming from a 6-quark. A kinematic fit is 
used to reconstruct the tt event topology. After back- 
ground subtraction, a Bayesian unfolding procedure is 
performed to correct for acceptance and detector ef- 
fects. The measured value of Ac is Ac = — 0.019 ± 
0.028 (stat.) ± 0.024 (syst.), consistent with the pre- 
diction from the MC@NLO Monte Carlo generator of 
Ac = 0.006±0.002. Measurements of Ac in two ranges 
of invariant mass of the top-antitop pair are also shown. 

Keywords Top physics • Charge asymmetry 

PACS 14.65.Ha 

1 Introduction 

The top quark is the heaviest elementary particle so far 
observed. With a mass close to the electroweak scale it 
may play a special role in physics beyond the Standard 
Model (SM) . Its pair production at hadron colliders al- 
lows a test of quantum chromodynamics (QCD) at high 
energies. 

This paper describes the measurement of the charge 
asymmetry Ac, defined as [1,2]: 
N{A\y\>0)-N{A\y\<0) 



A 



(1) 



7V(Z\|2/|>0)+iV(Z\|y|<0)' 
where A\y\ = \yt\~\yt\is the difference between the ab- 
solute values of the top and antitop rapidities {\yt\ and 
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\yt\) and N is the number of events with A\y\ positive 
or negative. 

Although tt production at hadron colliders is pre- 
dicted to be symmetric under the exchange of t and t at 
leading order, at next-to-leading order (NLO) the pro- 
cess qq — > ttg exhibits an asymmetry in the differential 
distributions of the top and antitop, due to interfer- 
ence between initial and final state gluon emission. The 
qq -^ ti process also possesses an asymmetry due to 
the interference between the Born and box diagrams. 
Similarly, the qg — >■ tiq process is asymmetric due to 
interference between amplitudes which have a relative 
sign difference under the exchange of t and i. The pro- 
duction of tt pairs by giuon-gluon fusion, gg — > tt, on 
the other hand, is symmetric. 

In pp collisions at the Tevatron, where top pairs are 
predominantly produced by quark- antiquark annihila- 
tion, perturbative QCD predicts that the top quark will 
be preferentially emitted in the direction of the incom- 
ing quark and the antitop in the direction of the in- 
coming antiquark [3]. Consequently, the charge asym- 
metry is measured as a forward-backward asymmetry, 
ApB- Recent measurements of ApB by the CDF and 
DO Collaborations [4-7] show a 2-3ct excess over the 
SM expectations enhancing interest in scrutinising the 
tt asymmetry. For tt invariant mass, mu, greater than 
450 CeV, the CDF experiment measures an asymme- 
try in the tt rest frame which is 3.4(t above the SM 
prediction [6]. Several new physics models have been 
proposed to explain the excess observed at CDF and 
DO [1,8-17]. Different models predict different asym- 
metries as a function of ttjjj [18]. 

In pp collisions at the LHC, the dominant mech- 
anism for tt production is expected to be the gluon- 
gluon fusion process, while tt production via qq or qg is 
small. Since the initial state is symmetric, the forward- 



backward asymmetry is no longer a useful observable. 
However, due to the asymmetry in the production via 
qq and qg, QCD predicts at the LHC a small excess 
of centrally produced antitop quarks while top quarks 
are produced, on average, at higher absolute rapidi- 
ties. This can be understood by the fact that for ti 
production via qq annihilation the valence quark car- 
ries, on average, a larger momentum fraction than the 
anti-quark from the sea. With top quarks preferentially 
emitted in the direction of the initial quarks in the tt 
rest frame, the boost into the laboratory frame drives 
the top mainly in the forward or backward directions, 
while antitops are preferentially retained in the central 
region. If new physics is responsible for the Tevatron 
ApB excess, the charge asymmetry measured at the 
LHC is a natural place to look for it. 

In this paper, the measurement of the charge asym- 
metry Ac is performed using candidate ti events se- 
lected in the lepton+jets channel. In this channel, the 
SM decay of the ti pair to W^hW^h results in a sin- 
gle electron or muon from one of the W boson decays 
and four jets, two from the second W boson decay and 
two from the h- and 6-quarks. To allow comparisons 
with theory calculations, the measured A\y\ distribu- 
tion is unfolded to account for acceptance and detector 
effects. An inclusive measurement, and measurements 
of Ac in two ranges of ii invariant mass, are presented. 
An inclusive measurement of this asymmetry with an 
equivalent observable has been recently reported by the 
CMS collaboration [19]. 



2 The ATLAS detector 

The ATLAS detector [20] at the LHC covers nearly the 
entire solid angle^ around the collision point. It consists 
of an inner tracking detector surrounded by a thin su- 
perconducting solenoid, electromagnetic and hadronic 
calorimeters, and an external muon spectrometer in- 
corporating three large superconducting toroid magnet 
assemblies. 

The inner-detector system is immersed in a 2 T axial 
magnetic field and provides charged particle tracking in 
the range \'q\ < 2.5. The high-granularity silicon pixel 
detector covers the vertex region and provides typically 



^ ATLAS uses a right-handed coordinate system with its ori- 
gin at the nominal interaction point (IP) in the centre of 
the detector and the z-axis along the beam pipe. The a;-axis 
points from the IP to the centre of the LHC ring, and the 
y axis points upward. Cylindrical coordinates (r, </>) are used 
in the transverse plane, <f) being the azimuthal angle around 
the beam pipe. The pseudorapidity is defined in terms of the 
polar angle 6 as -q = — lntan(6/2). Transverse momentum 
and energy are defined as px = psinS and Et = EsinO, 
respectively. 



three measurements per track, followed by the silicon 
microstrip tracker (SCT) which provides four measure- 
ments from eight strip layers. These silicon detectors 
are complemented by the transition radiation tracker 
(TRT), which enables extended track reconstruction up 
to 1 77 1 = 2.0. In giving typically more than 30 straw-tube 
measurements per track, the TRT improves the inner 
detector momentum resolution, and also provides elec- 
tron identification information. 

The calorimeter system covers the pseudorapidity 
range \r]\ < 4.9. Within the region \r]\ < 3.2, electro- 
magnetic calorimetry is provided by barrel and end- 
cap lead-liquid argon (LAr) electromagnetic calorime- 
ters, with an additional thin LAr presampler covering 
1 77 1 < 1.8 to correct for energy loss in material upstream 
of the calorimeters. Hadronic calorimetry is provided by 
the stecl/scintillating-tile calorimeter, segmented into 
three barrel structures within \ri\ < 1.7, and two cop- 
pcr/LAr hadronic endcap calorimeters. The solid an- 
gle coverage is completed with forward coppcr/LAr and 
tungsten/LAr calorimeter modules optimised for elec- 
tromagnetic and hadronic measurements respectively. 

The muon spectrometer comprises separate trigger 
and high-precision tracking chambers measuring the de- 
flection of muons in a magnetic field with a bending in- 
tegral from 2 to 8 Tm in the central region, generated 
by three superconducting air-core toroids. The preci- 
sion chamber system covers the region \r]\ < 2.7 with 
three layers of monitored drift tubes, complemented by 
cathode strip chambers in the forward region, where the 
background is highest. The muon trigger system covers 
the range \ri\ < 2.4 with resistive plate chambers in the 
barrel, and thin gap chambers in the endcap regions. 

A three-level trigger system is used to select in- 
teresting events. The level-1 trigger is implemented in 
hardware and uses a subset of detector information 
to reduce the event rate to a design value of at most 
75 kHz. This is followed by two software-based trigger 
levels, level-2 and the event filter, which together reduce 
the event rate to about 300 Hz. 

3 Data and Monte Carlo samples 

Data from LHC pp collisions collected by the ATLAS 
detector between March and June 2011 are used in the 
analysis, corresponding to an integrated luminosity of 
1.04 fb-i. 

Simulated top pair events are generated using the 
MCONLO [21] Monte Carlo (MC) generator with the 
NLO parton density function (PDF) set CTEQ6.6 [22]. 
Parton showering and the underlying event are mod- 
elled using HERWIG [23] and JIMMY [24], respec- 
tively. This ti sample is normalised to a cross section 



of 165 pb, obtained with the latest theoretical com- 
putation, which approximates the next-to-next-to lead- 
ing order prediction [25]. Single top events are also 
generated using MC@NLO while the production of 
W/Z bosons in association with jets is simulated us- 
ing the ALPGEN generator [26] interfaced to HER- 
WIG and JIMMY with CTEQ6.1 [27]. Diboson events 
{WW, WZ, ZZ) are generated using HERWIG with 
MRST20071omod [28]. 

All Monte Carlo simulation samples are generated 
with multiple pp interactions per bunch crossing (pile- 
up) . These simulated events are re- weighted so that the 
distribution of the number of interactions per crossing 
in simulation matches that in the data. The samples 
are then processed through the GEANT4 [29] simula- 
tion [30] of the ATLAS detector and the standard re- 
construction software. 



4 Event selection 

4.1 Physics object selection 

Reconstructing top quark pair events in the detector 
requires electrons, muons, jets and missing momentum 
to be simultaneously measured. Electron candidates are 
defined as energy deposits in the electromagnetic calorime- 
ter associated with a well-measured track. Identification 
criteria based on shower shape variables, track quality, 
and information from the transition radiation tracker 
are applied to electron candidates [31]. All candidates 
are required to have pt > 25 GeV and |??ciustor| < 2.47, 
where Tydustcr is the pscudorapidity of the electromag- 
netic calorimeter cluster associated with the electron. 
Candidates in the calorimeter transition region 1.37 < 
|??ciuster| < 1.52 are excluded. 

Muon candidates are reconstructed from track seg- 
ments in different layers of the muon chambers. These 
segments are combined starting from the outermost layer, 
with a procedure that takes material effects into ac- 
count, and matched with tracks found in the inner de- 
tector. The candidates are then refitted using the com- 
plete track information from both detector systems, and 
are required to satisfy pt > 20 GeV and I77I < 2.5. 

Jets are reconstructed with the anti-kt algorithm, 
with a distance parameter of 0.4 [32], starting from 
clusters of energy in adjacent calorimeter cells at the 
electromagnetic (EM) scale. The jet energy is corrected 
to the hadronic scale using p^- and 77-dependent cor- 
rection factors obtained from simulation and validated 
with data [33]. Jet quality criteria are applied to iden- 
tify jets not associated to in-time real energy deposits in 
the calorimeters caused by various sources (calorimeter 



noise, non-collision beam-related background, cosmic- 
ray induced showers). 

The missing transverse momentum {E^^^'^) is recon- 
structed from clusters of energy calibrated at the EM 
scale and corrected according to the energy scale of 
the associated physics object [34]. Contributions from 
muons are included using their momentum measured 
from the tracking and muon spectrometer systems. The 
remaining clusters not associated with the high px ob- 
jects are also included in the missing transverse mo- 
mentum. 

Muons within AR = 0.4 of a jet axis^ and with 
Pt > 20 GeV are removed in order to reduce the con- 
tamination caused by muons from hadron decays. Sub- 
sequently, jets within AR = 0.2 of an electron candidate 
are removed to avoid double counting electrons as jets. 

Isolation criteria are applied to both electron and 
muon candidates to reduce the backgrounds from hadrons 
mimicking lepton signatures and backgrounds from heavy 
flavour decays inside jets. For electrons, the total energy 
in a cone of AR = 0.2 around the electron candidate 
must not exceed 3.5 GeV, after correcting for energy de- 
posits from pile-up and for the energy associated with 
the electron. For muons, the sum of track transverse 
momenta for all tracks with px > 1 GeV and the to- 
tal energy deposited in a cone of AR = 0.3 around the 
muon are both required to be less than 4 GeV ignoring 
the contribution of the muon px- 

Reconstructing top quark pair events is facilitated 
by the ability to tag jets from the hadronisation of b- 
quarks. For this purpose, two 6-tagging algorithms are 
used and their results are combined to extract a tagging 
decision for each jet. One 6-tagger exploits the topol- 
ogy of b- and c-hadron weak decays inside the jet. A 
Kalman filter [35] is used to find a common line on 
which the primary vertex and the b- and c-hadron decay 
vertices lie, as well as their position on this line, giv- 
ing an approximate flight path for the b- and c-hadrons. 
The discrimination between 6-, c- and light quark jets is 
based on a likelihood using the masses, momenta, flight- 
length significances, and track multiplicities of the re- 
constructed vertices as inputs. To further increase the 
flavour discrimination power, a second 6-tagger is run 
which docs not attempt to directly reconstruct decay 
vertices. Instead, this second tagger uses the transverse 
and the longitudinal impact parameter significances of 
each track within the jet to determine a likelihood that 
the jet originates from a 6-quark. The results of both 
taggers are combined using a neural network to deter- 
mine a single discriminant variable which is used to 
make tagging decisions. The combined tagger operating 



^ AR = \J Af^"^ -\- At)^, where A(p and Ar] are the separation 
in azimuthal angle and pseudorapidity, respectively. 



point chosen for the present analysis corresponds to a 
70% tagging cfSciency for 5-jets in simulated tt events 
while light flavor jets are suppressed by approximately 
a factor of 100. 



4.2 Selection of tt candidates 

The tt final state in the lepton+jets channel is charac- 
terised by an isolated lepton (electron or muon) with 
relatively high px, missing transverse momentum aris- 
ing from the neutrino from the leptonic W decay, two 
6-quark jets and two light quark jets from the hadronic 
W decay. To select events with this topology, the appro- 
priate single-electron or single-muon trigger is required 
to have fired (with thresholds at 20 and 18 GeV respec- 
tively). The events are also required to contain one and 
only one reconstructed lepton with pT > 25 GeV for 
electrons and pt > 20 GeV for muons. To reject multi- 
jet background in the muon channel, E!f^^^ > 20 GeV 
and E^^^^ + mT{W) > 60 GeV are required-^. In the 
electron channel more stringent cuts on E!f^^^ and tot(W^) 
are required because of the higher level of multijet back- 
ground, i.e. Sf?'"" > 35 GeV and mT{W) > 25 GeV. 
Events are required to have at least four jets with px > 
25 GeV and |?/| < 2.5. These requirements define the 
'pretag' selection. The 'tagged' selection requires, in ad- 
dition, at least one of the jets with p^ > 25 GeV and 
I?;! < 2.5 to be 6-tagged. 



5 Background determination 

5.1 Multijet background 

The method used for evaluating the multijet background 
with fake leptons^ in both the electron and muon chan- 
nels is the so-called 'Matrix Method'. This relies on 
defining loose and tight lepton samples [36] and mea- 
suring the fractions of real (ej-cai) and fake (efakc) loose 
leptons that arc selected as tight leptons. The fraction 
eroai is measured using data control samples of Z boson 
decays to two leptons, while Cfakc is measured from data 
control regions defined separately for the electron and 
muon channels, where the contribution of fake leptons 
is dominant. 

For the muon channel, the loose data sample is de- 
fined by removing the isolation requirements in the de- 



^Here rmiW) is the W^-boson transverse mass, defined as 
■')) where the measured B™'^^ vector 



^ 



provides the neutrino information. 
*The term 'fake' leptons here refers to hadrons mimicking 
lepton signatures and to leptons arising from heavy hadron 
decays, whereas 'real' leptons come from W and Z decays. 



fault muon selection. The fake lepton efficiencies are de- 
termined using a low ttit control region tot < 20 GeV 
with an additional cut E'^""" -K tot < 60 GeV. The ef- 
ficiencies for signal and fake leptons are parameterised 
as a function of muon |r/| and pT in order to account for 
the variation of the muon detector acceptance and the 
profile of hadronic activity in the detector that affects 
the muon isolation. 

For the multijet background estimate in the electron 
channel, the loose data sample is defined by considering 
events with electrons passing looser identification crite- 
ria. The electron isolation requirement is also modified: 
the total energy in a cone of AR = 0.2 around the elec- 
tron is required to be smaller than 6 GeV (instead of 3.5 
GeV) , after correcting for energy deposits from pile- up 
interactions and for the energy associated with the elec- 
tron. The fake lepton efficiencies are determined using 
a low £:™'"" control region (5 GeV < E!^''''' < 20 GeV). 

In both channels contributions from VF-l-jets and 
Z-|-jets backgrounds in the control region, estimated 
using Monte Carlo simulation, are subtracted. 



5.2 M^-t-jets background estimation 

At the LHC the rate of iy"'"-|-jets is larger than that 
of iy~-|-jets because there are more valence u quarks 
than d quarks in the proton. Theoretically, the ratio of 
VF^-l-jets and M^~-|-jets cross sections is predicted much 
more precisely than the total M^-|-jets cross section [37, 
38]. This asymmetry is exploited here to measure the 
total ly-fjcts background from the data. 

Since, to a good approximation, processes other than 
M^-l-jets give equal numbers of positively and negatively 
charged leptons, the formula 

/ rAic + 1 ' 



N, 



w+ 



+ Nw- 



{D+-D-1 



(2) 



\rMC - 1, 

can be used to estimate the total number of W events in 
the selected sample. Here D^{D^) are the total num- 
bers of events in data passing the selection cuts de- 
scribed in Section 4.2 (apart from the 6-tagging require- 
ment) with positively (negatively) charged leptons, and 



tmc 



_ jv(pp-»ty+) 

N{pp^W-) 



is evaluated from Monte Carlo sim- 



ulation, using the same event selection. 

The ratio tmc is found to be 1.56 ±0.06 in the elec- 
tron channel and 1.65 ±0.08 in the muon channel. The 
dominant uncertainties on tmc originate from those of 
the parton distribution functions, the jet energy scale, 
and the heavy flavour fractions in ly+jets events (frac- 
tions of ly+jets events containing 66 pairs, cc pairs and 
c quarks). 

Since the theoretical prediction for heavy flavour 
fractions in VK±jets suffers from large uncertainties, a 



data-driven approach was developed to constrain these 
fractions with some inputs from MC simulation. In this 
approach samples with a lower jet multiplicity, obtained 
from the selection described in Section 4.2, but requir- 
ing precisely one or two jets instead of four or more 
jets, are analysed. The numbers VFj^retag' ^ulggod' o^ 
W+i jet events in these samples (where i = 1,2), before 
and after applying the 6-tagging requirement, are com- 
puted by subtracting the small contributions of other 
Standard Model processes - electroweak {WW, WZ, 
ZZ and Z-|-jets) and top {ti and single top) using pre- 
dictions from the simulation, and by subtracting the 
multijet background as described in Section 5.1. 

A system of two equations, expressing the number of 
T^ -I- 1 jet events and W + 2 jets events before and after 
6-tagging, can be written with six independent flavour 
fractions as the unknowns, corresponding to fractions 
of Wbb+iets, Tycc-l-jets, and T4^c-|-jets events in the one 
and two jet bins. The simulation prediction for the ra- 
tio of the heavy flavour fractions between the one and 
two jet bins is used to relate the heavy flavour fractions 
in the two bins, reducing the number of independent 
fractions to three. Finally, the ratio of the fractions of 
Wcc-|-jets and W^66-|-jets events in the two-jet bin is 
taken to be fixed to the value obtained from simulated 
events in order to obtain two equations for two indepen- 
dent fractions. Based on this measurement, the heavy 
flavour fractions in simulated M^-fjets events are ad- 
justed by a scale factor 1.63 ± 0.76 for Wbb+iets and 
iycc-|-jets events and 1.11 ± 0.35 for M^c-fjets. When 
applied to the signal region, an additional 25% uncer- 
tainty on these fractions is added, corresponding to the 
uncertainty of the Monte Carlo prediction for the ratio 
of heavy flavour fractions in different jet multiplicities. 
The heavy flavour scale factors are applied to simulated 
W^-l-jets events throughout this paper, and the effect of 
their uncertainties on the value of rMC is evaluated. 

Using Equation (2), the total number of W^-l-jets 
events passing the event selection described in Section 4.2 
without requiring a ^-tagged jet, VF>4_protag, is evalu- 
ated to be 5400 ± 800 (stat. -I- syst.) in the electron 
channel and 8600 ± 1200 (stat. -|- syst.) in the muon 
channel. 

The number of V7-|-jets events passing the selection 
with at least one 6-tagged jet is subsequently evaluated 
as [36] 



W^>4, tagged = W> 



4,prctag 



/2, tagged ■ k2^>4 



(3) 



Here /s.tagged = W'a^figged/W'i^pretag is the fraction of 
W+2 jets events passing the requirement of having at 
least one 6-tagged jet, and fc2^>4 = />i;tagged//2'laggod 
is the ratio of the fractions of simulated VF-|-jets events 
passing the requirement of at least one ^-tagged jet. 



for at least four and two jets, respectively. The value 
of /2, tagged is fouud to be 0.065 ± 0.005 in the elec- 
tron and 0.069 ± 0.005 in the muon channel, where the 
uncertainties include statistical and systematic contri- 
butions. The ratio k2^>4 is found to be 2.52 ± 0.36 in 
the electron channel and 2.35 ± 0.34 in the muon chan- 
nel. The uncertainties include both systematic contri- 
butions and contributions arising from the limited num- 
ber of simulated events. The total number of Ty-|-jets 
events passing the selection with a 6-tagged jet, W^>4,taggc 
is evaluated to be 880 ± 200 (stat. + syst.) in the elec- 
tron channel and 1390 ± 310 (stat. -t- syst.) in the muon 
channel. 



5.3 Other backgrounds 

The numbers of background events coming from single 
top production, Z-t-jets and diboson events are evalu- 
ated using Monte Carlo simulation normalized to the 
relevant NNLO cross sections for Z+jets events and 
NLO for diboson events. 



5.4 Event yield 

The final numbers of expected and observed data events 
in both channels after the full event selection are listed 
in Table 1 . The number of events in the electron chan- 
nel is significantly lower than in the muon channel due 
to the higher lepton px requirement and the more strin- 
gent missing momentum requirement, which are neces- 
sary to reduce the contribution from the multijet back- 
ground. The overall agreement between expectation and 
data is good. 



6 Reconstruction of the tt final state 

To measure the charge asymmetry in top pair events, 
the full tt system is reconstructed. For this purpose, a 
kinematic fit is used that assesses the compatibility of 
the observed event with the decays of a top-antitop pair 
based on a likelihood approach. 

The likelihood takes as inputs the measured ener- 
gies, pseudorapidities and azimuthal angles of four jets, 
the measured energy of the lepton, and the missing 
transverse momentum. If there are more than four jets 
in the event satisfying px > 25 GeV and \ri\ < 2.5, all 
subsets of four jets from the five jets in the event with 
highest pt are considered. 

The likelihood is computed as 

L = B{Ep^i, Ep^2\mw , Tw) ■ B{Eiep,E^\mw,rw) ■ 



Channel 



/i + jets pretag I fi + jets tagged 



e + jets pretag 



jets tagged 



tt 


7200 


± 


600 


6300 


± 


500 


4800 


± 


400 


4260 


± 


350 


W+jets 


8600 


± 


1200 


1390 


± 


310 


5400 


± 


800 


880 


± 


200 


Single top 


460 


± 


40 


366 


± 


32 


320 


± 


28 


256 


± 


22 


Z+jets 


940 


± 


330 


134 


± 


47 


760 


± 


270 


110 


± 


40 


Diboson 


134 


± 


7 


22 


± 


2 


80 


± 


5 


13 


± 


1 


Multijets 


1500 


± 


800 


500 


± 


500 


900 


± 


500 


250 


± 


250 



Total background 


11700 ± 


1400 


2400 ± 


600 


7500 ± 


900 


1500 ± 


320 


Signal + background 


18900 ± 


1600 


8800 ± 


800 


12000 ± 


1000 


5800 ± 


500 1 




Observed 


19639 




9124 




12096 




5829 





Table 1 Numbers of events observed in data and expected from ti signal events and various background processes for the 
pretag and tagged samples defined in Section 4.2. The experimentally determined uncertainties quoted for VK+jets and multijet 
backgrounds include systematic uncertainties on the normalisation. The quoted uncertainties on the other backgrounds are 
those from theory, taken to be 8% for tt and single top, 34% for Z+jets and 5% for diboson backgrounds. The numbers 
correspond to an integrated luminosity of 1.04 fb~^ in both electron and muon channels. 



B{Ep^i,Ep^2, Ep^3\mt,rt) ■ B{Eicp, E^, Ep^4\mt,rt) 

W{Ef''\p,,,) ■ yV{Ef''\py,,) ■ >V(i^lep|^lep) • 
4 4 

Y[WiEjct,i\Eps) • J|P(taggcd | parton flavour), 

i=l i=l 

where: 



Distributions of the invariant mass and transverse 
momentum of the reconstructed top-antitop pair are 
shown in Fig. 1. 



Symbols B represent Breit-Wigncr functions, eval- 
uated using invariant masses of sums of appropri- 
ate parton and lepton four-vectors. The pole masses 
of the W boson and the top quark are fixed to 
mw = 80.4 GeV and nit = 172.5 GeV, respectively. 
Their widths are taken to be /V = 2.1 GeV and 
A = 1.5 GeV. 

Symbols W represent the transfer functions asso- 
ciating the reconstructed quantities {X) to quarks 
and leptons produced in the hard scattering {X). 
Ep^i are the energies of partons associated to jets Si = 2^ RijTj 



7 Unfolding 

The measured distributions of top and anti-top rapidi- 
ties are distorted by detector effects and an event se- 
lection bias. To correct for these distortions the experi- 
mental distributions are unfolded to the four-vectors of 
the top quarks before decay. 

The relation between a true distribution Tj (assum- 
ing, for simplicity, that there is only one observable of 
interest) and the reconstructed distribution Si after de- 
tector simulation and event selection can be written: 



(4) 



with measured energies -Ejet.i- These transfer func- 
tions are derived from Monte Carlo simulation. 
— P(tagged I parton flavour) is the ^-tagging probabil- 
ity or rejection efficiency, depending on the parton 
flavour, as obtained from Monte Carlo simulation. 

The likelihood is maximised with respect to the en- 
ergies of the partons, the energy of the charged lepton, 
and the components of the neutrino three-momentum. 
The assignment of jets to partons which gives the high- 
est likelihood value is selected. Finally, the sign of the 
charge of the top quark (or anti-quark) decaying into 
the lepton is determined from the lepton charge. 

The overall efficiency for the reconstruction of the 
correct event topology is found to be 74% in Monte 
Carlo simulated ii events. Only those events where four 
jets and a lepton are matched to partonic particles are 
considered for the efficiency computation. 



where Rij is the response matrix defined as the proba- 
bility to observe an event in bin i when it is expected 
in bin j. 

The true distribution Tj can be obtained from the 
observed distribution 5"^ by inverting the response ma- 
trix. The unfolding problem can similarly be formulated 
for the case of multiple observables. In this analysis, 
Bayes' theorem is applied iteratively in order to per- 
form the unfolding [39]. 

The unfolding is performed using response matrices 
which account for both detector response and accep- 
tance effects. The response matrices are calculated us- 
ing Monte Carlo events generated with MC@NLO. The 
unfolding is done separately, after background subtrac- 
tion, for the inclusive measured distribution of zl|j/| (a 
one-dimensional unfolding problem), and the measured 
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Fig. 1 Expected and observed distributions for the invariant mass (plots (a) and (b)) and transverse momentum (plots (c) 
and (d)) of the reconstructed ti system. The left hand panels show distributions in the electron channel, while the right 
hand panels show distributions in the muon channel. The data are compared to the sum of the ti signal contribution and 
backgrounds. The background contributions from VK+jets and multijet production have been estimated from data, while the 
other backgrounds are estimated from simulation. The uncertainty on the combined signal and background estimate includes 
systematic contributions. Overflows are shown in the highest bin of each histogram. 



distribution A\y\ as a function of the reconstructed top- 
antitop invariant mass ijif-i (a two-dimensional unfold- 
ing problem). 



Two bins are used for m^j in the two-dimensional 
unfolding of zi|y| versus rriu, separated at mti = 450 GeV. 
The choice of this rritt value is motivated by the ob- 
served CDF forward-backward asymmetry [6] and by 
separating the data sample into two bins with roughly 
equal number of events. 



An additional cut on the value of the likelihood for 
the tt candidate is required in the two-dimensional un- 
folding, since a large fraction of simulated events with 
a badly reconstructed ?7i,(t are found to have a low like- 
lihood value. 

The response matrix (including both detector and 
acceptance effects) for the inclusive Ac measurement 
is shown in Fig. 2. Six bins in A\y\, in the range — 3 < 
A\y\ < 3, are used in the response matrix, with the 
outermost bins broader than the inner bins in order 



to avoid the occurrence of bins with no entries in the 
measured distributions. Only a very small fraction of 
simulated tt events are found to have |^|?/|| > 3, and 
hence such events have a negligible influence on the 
results. 

The unfolding procedure is applied to the observed 
A\y\ distribution in data, after subtracting background 
contributions. When performing the background sub- 
traction, the shape of the multijet background is ob- 
tained by applying the Matrix Method (described in 
Section 5.1) in bins of A\y\. The shape of all remain- 
ing backgrounds is taken from Monte Carlo simulation. 
The value of Ac after unfolding is obtained by counting 
the numbers of events with A\y\ > and A\y\ < in 
the unfolded A\y\ distribution. 



8 Systematic uncertainties 

Several sources of systematic uncertainties are taken 
into account in this analysis. These arc categorised into 
the detector modelling, the modelling of signal and back- 
ground processes and the unfolding method. 



8.1 Detector modelling 

Small mis-modellings of muon or electron trigger, recon- 
struction and selection efficiencies in simulation are cor- 
rected for by scale factors derived from measurements of 
the efficiency in data. Z — > fxfi or Z — >■ ee and W —^ ev 
decays are used to obtain scale factors as functions of 
the lepton kinematics. The uncertainties are evaluated 
by varying the lepton and signal selections and from the 
uncertainty in the evaluation of the backgrounds. Sys- 
tematic uncertainties at the level of 1% are found for 
both cases. The same processes are used to measure the 
lepton momentum scale and resolution. Scale factors, 
with uncertainties at the level of (1-1.5)%, are derived 
to match the simulation to observed distributions. A 
systematic uncertainty for charge mis-identification of 
leptons is assigned which is negligible for muons and 
ranges from 0.2% to 3% for electrons depending on \rj\. 
The jet energy scale is derived using information 
from test-beam data, collision data and simulation. Its 
uncertainty varies between 2.5% and 8% in the central 
region, depending on jet px and 77 [33]. This includes 
uncertainties in the flavour composition of the sample 
and mis-measurements due to the effect of nearby jets. 
Pile-up gives additional uncertainties of up to 5% (7%) 
in the central (forward) region. An extra uncertainty of 
0.8% to 2.5%, depending on jet px, is assigned to jets 
arising from the fragmentation of 6-quarks, due to dif- 
ferences between light and gluon jets as opposed to jets 



containing 6-hadrons. The jet energy resolution and re- 
construction efficiency are measured in data using tech- 
niques described in Rcfs. [33,40], and their uncertainties 
are found to be 10% and (1-2)%, respectively. 

The 6-tagging efficiencies and mis-tag rates are mea- 
sured in data. Jet px dependent scale factors, applied to 
simulations to match the efficiencies measured in data, 
have uncertainties which range from 9% to 15% and 
11% to 22%, respectively. A systematic uncertainty is 
assigned for a potential difference of up to 5% between 
the 6-tagging efficiency for 5-jets and that of ti-jets. The 
uncertainty on the measured luminosity is 3.7% [41,42]. 

Due to a hardware failure, later repaired, one small 
region of the liquid argon calorimeter could not be read 
out in a subset of the data corresponding to 84% of 
the total integrated luminosity. Data events in which 
an electron or jet with px > 20 GeV is close to the 
affected calorimeter region are rejected for the rele- 
vant part of the dataset. Monte Carlo simulated events 
with electrons or jets of px > 20 GeV close to the af- 
fected region are rejected with a probability equal to 
the fraction of the integrated luminosity of data for 
which the calorimeter hardware problem was present. 
A systematic uncertainty is evaluated by varying the 
Px-threshold in data of the electrons and jets near the 
affected region by ±4 GeV, corresponding to the un- 
certainty in the energy lost by objects in the affected 
region. 



8.2 Signal and background modelling 

The systematic uncertainty in the modelling of the sig- 
nal process is assessed by simulations based on dif- 
ferent Monte Carlo generators. Sources of systematic 
uncertainty considered here are the choice of gener- 
ator and parton shower model, the choice of parton 
density functions, the assumed top quark mass and 
the choice of parameters which control the amount of 
initial and final state radiation. Predictions from the 
MCONLO and POWHEG [43,44] generators arc com- 
pared. The parton showering is tested by comparing 
two POWHEG samples interfaced to HERWIG and 
PYTHIA, respectively. The amount of initial and fi- 
nal state radiation is varied by modifying parameters 
in ACERMC [45] interfaced to PYTHIA according to 
Ref. [46]. The parameters are varied in a range com- 
parable to those used in the Perugia Soft/Hard tune 
variations [47] . The impact of the choice of parton den- 
sity functions is studied using the procedure described 
in Ref. [48]. MC@NLO samples are generated assum- 
ing different top quark masses and their predictions are 
compared. The observed differences in the results are 
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Fig. 2 Correlations between tlie true and reconstructed values of A\y\ encoded in the unfolding response matrix for the 
electron (left) and muon (right) channels. The value of an entry in the matrix is proportional to the area of the corresponding 
box. 



scaled to variations of ±0.9 GeV according to the un- 
certainty on the measured value [49] . 

As described in Section 5, background processes are 
either modelled by simulation or estimated in auxiliary 
measurements. The uncertainty in the estimate of the 
multijet background is evaluated by considering mod- 
ified definitions of the loose data sample, taking into 
account the statistical uncertainty in measuments of 
Croaijefake described in Section 5.1 as well as the un- 
certainties in the normalisations of the M^-|-jets and 
Z+icts backgrounds which arc subtracted in the con- 
trol region. The total uncertainty is estimated to be 
100%. The normalisation of M^+jets processes is evalu- 
ated from auxiliary measurements using the asymmet- 
ric production of positively and negatively charged W 
bosons in VF-|-jets events. The uncertainty is estimated 
to be 21% and 23% in the four jet bin, for the elec- 
tron and muon channels respectively. This uncertainty 
was estimated by evaluating the effect on both Vmc 
and fc2-i.>4 from the JES uncertainty and different PDF 
and generator choices. Systematic uncertainties on the 
shape of M^-fjets distributions are assigned based on 
differences in simulated events generated with differ- 
ent simulation parameters. Scaling factors correcting 
the fraction of heavy flavour contributions in simulated 
W^-|-jets samples are estimated in auxiliary measure- 
ments, as described in Section 5.2. The systematic un- 
certainties are found by changing the normalisations 
of the non-W processes within their uncertainties when 
computing V^i^^otag' ^i^llged' ^^ '^^11 ^^ taking into ac- 
count the impact of uncertainties in 6-tagging efficien- 



cies. The total uncertainties are 47% for H^&6-|-jets and 
VFcc-|-jets contributions and 32% for VFc-|-jets contribu- 
tions. The normalisation of Z+^et events is estimated 
using Berends-Giele-scaling [50] . The uncertainty in the 
normalisation is 48% in the four jet bin and increases 
with the jet multiplicity. A systematic uncertainty in 
the shape is accounted for by comparing simulated sam- 
ples generated with ALPGEN and SHERPA [51]. The 
uncertainty on the normalisation of the small back- 
ground contributions from single top and diboson pro- 
duction is estimated to be about 10% (depending on 
the channel) and 5%, respectively. 

Limited Monte Carlo sample sizes give rise to a sys- 
tematic uncertainty in the response matrix. This is ac- 
counted for by independently varying the bins of the 
response matrix according to Poisson distributions. 



8.3 Uncertainties from unfolding 

Closure tests are performed in order to check the valid- 
ity of the unfolding procedure. Reweighted tt samples 
with different amounts of asymmetry are considered. 
Pseudoexperiments are performed, varying the entries 
in histograms of the reconstructed distribution, to con- 
firm that the response of the unfolding is linear in the 
true value of Ac and that the true value of Ac is re- 
covered on average. A total of 40 iterations are used in 
both channels for the inclusive Ac measurement. For 
the measurement of Ac as a function of r/i^f, 80 itera- 
tions are used. The number of iterations is chosen by 
ensuring that the unfolding procedure has converged in 
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the sense that the absolute change in the unfolded value 
of Ac after performing an extra iteration is less than 
0.001. It is found that the unfolded values of Ac from 
all pscudoexperimcnts and the data converge before the 
chosen numbers of iterations. The potential bias aris- 
ing from the choice of convergence criterion is taken 
into account by adding an additional systematic un- 
certainty corresponding to the change in the unfolded 
value of Ac obtained by further increasing the number 
of iterations to very large values (10^). 

Pull distributions are constructed from pseudocx- 
periments and a relative shift of between 0% and 10% 
is found in the unfolded value of Ac with respect to 
the true value. An extra systematic uncertainty is as- 
signed to the unfolded value of Ac obtained from data, 
corresponding to this shift. 

In pseudoexperimcnts, a small bias is observed in 
the unfolded distributions corresponding to a relative 
difference of a few percent between the unfolded result 
and true value in each bin. An additional relative un- 
certainty of (2-5)% is applied to all bins of the unfolded 
distributions, corresponding to the largest relative bin 
deviation observed in pseudoexperimcnts. 

The statistical uncertainty in the unfolded measure- 
ment was computed using pscudoexperimcnts, propa- 
gating the uncertainties from the measured distribution 
using the statistical correlation matrix. 



8.4 Impact of systematic uncertainties 

The impact of the systematic uncertainties is evalu- 
ated by modifying the subtracted background before 
unfolding and by modifying the response matrix used 
for unfolding when relevant. In particular the detec- 
tor modelling systematic uncertainties are evaluated by 
shifting the estimated background as well as modifying 
the response matrix. Signal modelling uncertainties are 
computed by replacing the response matrix, and back- 
ground modelling uncertainties by modifying the esti- 
mated background. 

Table 2 summarises the sources of systematic un- 
certainties for the inclusive measurement of the charge 
asymmetry, and their impact on the measured asym- 
metry, after unfolding. The systematics for the two rriff 
bins are determined in a similar fashion. The evaluation 
of some systematic uncertainties is limited by the finite 
size of the Monte Carlo samples. In these cases, the 
larger of the electron and muon channel uncertainties 
is used for the uncertainty on the combined result. The 
resulting combined systematic uncertainties are ±0.028 
in the electron channel and ±0.024 in the muon chan- 
nel. 



9 Summary of Results 

The measured distributions of the top-antitop rapidity 
difference A\y\ = \yt\ ~ \yi\ before unfolding are shown 
in Fig. 3 for the electron and muon channel. Fig. 4 
shows the corresponding A\y\ distributions after un- 
folding. After unfolding, the bins of the measured dis- 
tribution have statistical and systematic correlations. 
Adjacent bins of the A\y\ distributions are found to 
be statistically anti-correlated with negative correlation 
coefficients of up to -0.6, whereas other correlations are 
small. 

The measured values of the top charge asymmetry 
before and after unfolding, defined by Eq. (1) in terms 
of A\y\, are summarised in Table 3. The analytic best 
linear unbiased estimator (BLUE) method [52, 53] is 
used to combine the measurement in the electron and 
muon channels after correction for detector resolution 
and acceptance. 

The measured asymmetries are: 

Ac = -0.019 ± 0.028 (stat.) ± 0.024 (syst.) 
for the integrated sample, and 

Ac = -0.052 ± 0.070 (stat.) ± 0.054 (syst.) 

for mti < 450 GeV, 
Ac = -0.008 ± 0.035 (stat.) ± 0.032 (syst.) 

for 77iti > 450 GeV. 

The measurement for the integrated sample can be com- 
pared with the result of the CMS Collaboration, Aq = 
-0.013 ±0.028 (stat)+o;°3^ (syst) [19]. Fig. 5 summa- 
rizes the measurements for the two m^t regions. These 
results are compatible with the prediction from the 
MC@NLO Monte Carlo generator of Ac = 0.006 ± 
0.002^, showing no evidence for an enhancement from 
physics beyond the Standard Model. 



10 Comparison of LHC and Tevatron Results 

The measurement of the charge asymmetry at the LHC 
is a test of the unexpectedly large forward-backward 
asymmetry observed at the Tevatron. However, because 
the LHC is a pp collider and the centre of mass energy is 
around three times larger, any relation between the two 
asymmetries is model-dependent. Here a comparison is 
made between the predicted values of the Tevatron and 



^The prediction of 0.0115 ± 0.0006 for the charge asymme- 
try found in Ref. [54] differs from the MC@NLO prediction 
of 0.006 ± 0.002, due to the former taking the LO prediction 
for the denominator in the definition (1) of Ac, and taking 
into account QED effects. The uncertainty on the MCQNLO 
prediction is obtained by considering variations in the renor- 
malisation and factorisation scales and different sets of PDFs. 
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Source of systematic uncertainty on Ac 


Electron channel 


Muon channel | 


Detector modelling | 


Jet energy scale 


0.012 


0.006 


Jet efficiency and resolution 


0.001 


0.007 


Muon efficiency and resolution 


<0.001 


0.001 


Electron efficiency and resolution 


0.003 


0.001 


b-tag scale factors 


0.004 


0.002 


Calorimeter readout 


0.001 


0.004 


Charge mis-fD 


<0.001 


<0.001 


b-tag charge 


0.001 


0.001 


Signal and background modelling | 


Parton shower/fragmentation 


0.010 


0.010 


Top mass 


0.007 


0.007 


tt modelling 


0.011 


0.011 


fSR and FSR 


0.010 


0.010 


PDF 


<0.001 


<0.001 


W+jets normalization and shape 


0.008 


0.005 


Z+jets normalization and shape 


0.005 


0.001 


Multijet background 


0.011 


0.001 


Single top 


<0.001 


<0.001 


Diboson 


<0.001 


<0.001 


MC Statistics 


0.006 


0.005 


Unfolding convergence 


0.005 


0.007 


Unfolding bias 


0.004 


<0.001 


Luminosity 


0.001 


0.001 



Total systematic uncertainty 



0.028 



0.024 



Table 2 List of sources of systematic uncertainties and their impact on the measured asymmetry in the electron and muon 
channel. In cases where asymmetric uncertainties were obtained, a symmetrisation of the uncertainties was performed by taking 
the average of the absolute deviations under systematic shifts from the nominal value. 
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Fig. 3 The measured A\y\ distribution before unfolding for the electron channel (left) and for the muon channel (right) after 
6-tagging is applied. Data (points) and Monte Carlo estimates (solid lines) are represented. The multijet background and the 
normalisation of the W+jets background are obtained as explained in Section 5. The uncertainty on the combined signal and 
background estimate includes both statistical and systematic contributions. 
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Fig. 4 The unfolded A\y\ distribution for tiie electron channel (left) and the muon channel (right) after f)-tagging, compared 
to the prediction from MC@NLO. The uncertainties on the measurement include both statistical and systematic contributions, 
which are shown separately. The inner part of the error bars corresponds to the statistical component of the uncertainty, while 
the outer part corresponds to the systematic component. The error bands on the MCQNLO prediction include uncertainties 
from parton distribution functions and renormalisation and factorisation scales. 



Asymmetry reconstructed detector and acceptance unfolded 


Ac (electron) -0.034 ± 0.019 (stat.) ± 0.010 (syst.) -0.047 ± 0.045 (stat.) ± 0.028 (syst.) 


Ac (muon) -0.010 ± 0.015 (stat.) ± 0.008 (syst.) -0.002 ± 0.036 (stat.) ± 0.024 (syst.) 


Combined -0.019 ± 0.028 (stat.) ± 0.024 (syst.) 



Table 3 The measured inclusive charge asymmetry values for the electron and muon channels after background substraction, 
before and after unfolding. 
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Fig. 5 Unfolded asymmetries in two regions of 171^^ compared to the prediction from MCQNLO. The error bands on the 
MC@NLO prediction include uncertainties from parton distribution functions and renormalisation and factorisation scales. 



LHC asymmetries for a few simpfe models beyond the 
SM. These are: (i) a flavour-changing Z' boson with 
right-handed couplings, exchanged in the t channel in 
uu — > ti [10]; (ii) a W boson, also with right-handed 



couplings, contributing in dd — >■ tt [11]; a heavy ax- 
igluon Qfj^ exchanged in the s channel [8,9]; (iv) a scalar 
doublet (j), with the same quantum numbers as the SM 
Higgs [55]; (v) a charge 4/3 scalar, colour-sextet (fi^) 
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or colour-triplet (w"*), contributing in the u channel to 
uu — )■ ti [12, 13]. In all these models, the parameter 
space is described by the mass M of the new particle 
(except for the axigluon which is assumed heavy, with 
M ^ 7 TeV) and a single coupling g. 

In order to find the correlated predictions for the 
forward-backward and charge asymmetries in each model, 
a comprehensive scan over the mass M and the cou- 
pling g is performed using the Protos generator [56], 
considering masses between 100 GeV and 10 TeV and 
the range of couplings for which the new physics con- 
tribution to the ti cross section at the Tevatron lies in 
the interval [-0.8,1.7] pb. This is a conservative require- 
ment which takes into account the different predictions 
for the SM cross section as well as the experimental 
measurement (see Rcf. [17] for details). 

In addition, a conservative upper limit on new physics 
contributions to o-jt for ttijj > 1 TeV is imposed. Fur- 
ther details can be found in Rcfs [17,55]. The coloured 
areas in Fig. 6 (a) represent the ranges of predicted val- 
ues for the inclusive Tevatron forward-backward asym- 
metry, ApB, and the inclusive LHC charge asymmetry, 
Ac, for the new physics models. The new physics con- 
tributions are computed using the tree-level SM ampli- 
tude plus the one(s) from the new particle(s). To a good 
approximation, the total asymmetries Aps , Ac are ob- 
tained from the former by summing the SM contribu- 
tion (at NLO in the lowest order). The horizontal lines 
correspond to the present ATLAS measurement and the 
measurement reported by the CMS Collaboration [19]. 
The vertical lines correspond to the asymmetry mea- 
surements at the Tevatron, ApB = 0.158 ± 0.075 [6] 
and Afb = 0.196 ±0.065 [7]. 

The ATLAS charge asymmetry measurement dis- 
favours models with a new flavour-changing Z' or W' 
vector boson proposed to explain the measured Teva- 
tron asymmetry. Minimal Z' models are also excluded 
by the non-observation of same-sign top quark produc- 
tion [57] . For the other new physics models the asymme- 
tries measured at the Tevatron are consistent with this 
measurement, within the experimental uncertainties. 

Fig. 6 (b) shows the allowed regions for the high- 
mass asymmetries (jritf > 450 GeV) at the Tevatron 
and the LHC for the six new physics models. The verti- 
cal lines represent the CDF measurement Aps = 0.475± 
0.114 [6], while the horizontal lines correspond to the 
present ATLAS measurement. In both panels of Fig. 6, 
the range of variation of SM predictions found in Refs [54, 
58,59] is indicated by a box. The predictions of the six 
new physics models are in tension with the CDF and 
ATLAS high- mass measurements considered together. 



11 Conclusion 

To summarize, the top quark charge asymmetry was 
measured in tt events with a single lepton (electron or 
muon), at least four jets and large missing transverse 
momentum using an integrated luminosity of 1.04 fb~^ 
recorded by the ATLAS experiment at a centre of mass 
energy of -^i = 7 TeV. The reconstruction of tt events 
was performed using a kinematic fit. The reconstructed 
inclusive distribution of A\y\ and the distribution as 
a function of m^f were unfolded after background sub- 
traction to obtain results that can be directly compared 
with theoretical computations. The results are compat- 
ible with the prediction from the MC@NLO Monte 
Carlo generator. These measurements disfavour models 
with a new flavour-changing Z' or W' vector boson that 
have been suggested to explain the measured Tevatron 
asymmetry. 
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Fig. 6 Measured FB asymmetries from the Tevatron and charge asymmetries from the LHC, compared to predictions from 
the SM as well as predictions incorporating various potential new physics contributions. The horizontal (vertical) bands and 
lines correspond to the ATLAS and CMS (CDF and DO) measurements. In (a) the inclusive values are presented and in (b) 
the ATLAS measurement for m^j > 450 GeV is compared to the CDF measurement. The MC predictions for the new physics 
models are from Refs. [17,55]. 



(Netherlands), PIC (Spain), ASGC (Taiwan), RAL (UK) 
and BNL (USA) and in the Tier-2 facihties worldwide. 



References 

10. 

1. S. Jung, A. Pierce, and J. D. Wells, Top quark 
asymmetry from a non-Abelian horizontal symmetry, 
Phys. Rev. D83 (2011) 114039, arXiv: 1103.4835 

[hep-ph]. 21. 

2. R. Diener, S. Godfrey, and T. A. W. Martin, Using 
Final State Pseudorapidities to Improve s-channel 
Resonance Observables at the LHC, Phys. Rev. D80 12. 
(2009) 075014, arXiv: 0909. 2022 [hep-ph]. 

3. J. H. Klihn and G. Rodrigo, Charge asymmetry in 
hadroproduction of heavy quarks, Phys. Rev. Lett. 81 

(1998) 49-52, arXiv:hep-ph/9802268. I3. 

4. DO Collaboration, V. M. Abazov et al.. Measurements 
of the Forward-Backward Charge Asymmetry in 
Top-Quark Pair Production, Phys. Rev. Lett. 100 
(2008) 142002, arXiv: 0712. 0851 [hep-ex]. 

5. CDF Collaboration, T. Aaltonen et al., 14. 
Forward-Backward Asymmetry in Top Quark Pair 
Production in pp collisions at \fs = 1.96 TeV , Phys. 

Rev. Lett. 101 (2008) 202001, arXiv: 0806. 2472 

[hep-ex] . 15. 

6. CDF Collaboration, T. Aaltonen et al., Evidence for a 
mass dependent forward-backward asymmetry in top 
quark pair production, Phys. Rev D83 (2011) 112003, 

arXiv: 1101.0034 [hep-ex]. 16. 

7. DO Collaboration, V. M. Abazov et al.. 
Forward-backward asymmetry in top quark- antiquark 
production, Phys. Rev. D84 (2011) 112005, 17. 
arXiv: 1107.4995 [hep-ex]. 



P. Ferrario and G. Rodrigo, Constraining heavy colored 

resonances from top-antitop quark events, Phys. Rev. 

D80 (2009) 051701, arXiv: 0906. 5541 [hep-ph]. 

P. H. Frampton, J. Shu, and K. Wang, Axigluon as 

Possible Explanation for pp — > ti Forward-Backward 

Asymmetry, Phys. Lett. B683 (2010) 294-297, 

arXiv:0911.2955 [hep-ph]. 

S. Jung, H. Murayama, A. Pierce, and J. D. Wells, Top 

quark forward-backward asymmetry from new t-channel 

physics, Phys. Rev. D81 (2010) 015004, 

arXiv:0907.4112 [hep-ph]. 

K. Cheung, W.-Y. Keung, and T.-C. Yuan, Top Quark 

Forward-Backward Asymmetry, Phys. Lett. B682 

(2009) , arXiv: 0908. 2589 [hep-ph]. 

J. Shu, T. M. Tait, and K. Wang, Explorations of the 

Top Quark Forward-Backward Asymmetry at the 

Tevatron, Phys. Rev. D81 (2010) 034012, 

arXiv:0911.3237 [hep-ph]. 

I. Dorsner, S. Fajfer, J. F. Kamenik, and N. Kosnik, 

Light colored scalars from grand unification and the 

forward-backward asymmetry in top quark pair 

production, Phys. Rev. D81 (2010) 055009, 

arXiv:0912.0972 [hep-ph]. 

B. Grinstein, A. L. Kagan, M. Trott, and J. Zupan, 

Forward-backward asymmetry in ti production from 

flavour symmetries, Phys. Rev. Lett. 107 (2011) 

012002, arXiv: 1102. 3374 [hep-ph]. 

Z. Ligeti, G. M. Tavares, and M. Schmaltz, Explaining 

the ti forward-backward asymmetry without dijet or 

flavor anomalies, JHEP 06 (2011) 109, 

arXiv: 1103.2757 [hep-ph]. 

J. A. Aguilar-Saavedra and M. Perez- Victoria, Probing 

the Tevatron ti asymmetry at LHC, JHEP 05 (2011) 

034, arXiv: 1103.2765 [hep-ph]. 

J. Aguilar-Saavedra and M. Perez- Victoria, 

Asymmetries in ti production: LHC versus Tevatron, 



15 



Phys.Rev. D84 (2011) 115013, arXiv: 1105.4606 
[hep-ph] . 

18. J. Aguilar-Saavedra and M. Perez- Victoria, Shaping the 
top asymmetry, arXiv: 1107.2120 [hep-ph]. 

19. CMS Collaboration, Measurem,ent of the charge 
asymmetry in top-quark pair production in 
proton-proton collisions at sqrt(s) = 7 TeV, Phys. Lett. 
B 709 (2012) 28, arXiv:1112.5100 [hep-ex]. 

20. ATLAS Collaboration, The ATLAS Experiment at the 
CERN Large Hadron Collider, JINST 3 (2008) S08003. 

21. S. Frixione and B.R. Webber, Matching NLO QCD 
computations and parton shower simulations, JHEP 
0206 (2002) 029, hep-ph/0204244. 

22. P. M. Nadolsky et al.. Implications of CTEQ global 
analysis for collider observables, Phys. Rev. D78 
(2008) 013004, arXiv: 0802. 0007 [hep-ph]. 

23. G. Corcella et al., HERWIG 6: An Event generator for 
hadron emission reactions with interfering gluons 
(including supersymmetric processes), JHEP 0101 
(2001) 010, hep-ph/0011363. 

24. J. M. Butterworth, J. R. Forshaw, and M. H. Seymour, 
Multiparton interactions in photoproduction at HERA, 
Z. Phys. C72 (1996) 637-646, arXiv:hep-ph/9601371. 

25. M. Aliev et al., HATHOR: HAdronic Top and Heavy 
quarks crOss section calculatoR, Comput. Phys. 
Commun. 182 (2011) 1034-1046, arXiv: 1007. 1327 
[hep-ph] . 

26. M.L. Mangano et al., ALPGEN, a generator for hard 
multiparton processes in hadronic collisions , JHEP 
0307 (2003) 001, hep-ex/0206293. 

27. Pumplin, J. and others. New generation of parton 
distributions with uncertainties from global QCD 
analysis, JHEP 07 (2002) 012, arXiv:hep-ph/0201195. 

28. A. Sherstnev and R. Thorne, Parton Distributions for 
LO Generators, Eur.Phys.J. C55 (2008) 553-575, 
arXiv:0711.2473 [hep-ph]. 

29. S. Agostinelli et al., Geant4 A Simulation Toolkit, 
Nucl. Instrum. Meth. A506 (2003) 250. 

30. ATLAS Collaboration, The ATLAS Simulation 
Infrastructure, European Physical Journal C - Particles 
and Fields 70 (2010) 823-874. 

31. ATLAS Collaboration, Electron performance 
measurements with the ATLAS detector using the 2010 
LHC proton-proton collision data, arXiv: 1110.3174 
[hep-ex] . 

32. M. Cacciari, G. P. Salam, and G. Soyez, The Anti-k(t) 
jet clustering algorithm, JHEP 0804 (2008) 063, 
arXiv:0802.1189 [hep-ph] . 

33. ATLAS Collaboration, Jet energy measurement with the 
ATLAS detector in proton-proton collisions at sqrt(s) 

= 7 Tel/, arXiv: 1112.6426 [hep-ex]. 

34. ATLAS Collaboration, Performance of Missing 
Transverse Momentum Reconstruction in 
Proton-Proton Collisions at 7 TeV with ATLAS, 
Eur.Phys.J. C72 (2012) 1844, arXiv: 1108.5602 
[hep-ex] . 

35. R. E. Kalnian, A new approach to linear filtering and 
prediction problems. Journal Of Basic Engineering 82 
(1960) no. Series D, 35-45. 

36. ATLAS Collaboration, Measurement of the top 
quark-pair production cross section with ATLAS in pp 
collisions at ^ = 7 TeV, Eur.Phys.J. C71 (2011) 
1577, arXiv: 1012. 1792 [hep-ex]. 

37. A. Martin, W. Stirling, R. Thorne, and G. Watt, 
Parton distributions for the LHC , The European 
Physical Journal C - Particles and Fields 63 (2009) 
189-285, arXiv: 0901. 0002 [hep-ph]. 



38. C.-H. Kom and W. Stirling, Charge asymmetry in W -t 
jets production at the LHC , The European Physical 
Journal C - Particles and Fields 69 (2010) 67-73, 
arXiv : 1004 . 3404 [hep-ph] . 

39. G. D'Agostini, A Multidimensional unfolding method 
based on Bayes' theorem, Nucl. Instrum. Meth. A362 
(1995) 487-498. 

40. ATLAS Collaboration, Jet energy resolution and 
selection efficiency relative to track jets from in-situ 
techniques with the ATLAS detector using 
proton-proton collisions at a center of mass energy ^/s 
= 7 TeV, ATLAS-CONF-2010-054 (2010) , 
http://cdsweb. cern.ch/record/1281311. 

41. ATLAS Collaboration, Luminosity Determination in pp 
Collisions at sqrt(s)=7 TeV Using the ATLAS Detector 
at the LHC, Eur.Phys.J. C71 (2011) 1630, 

arXiv: 1101.2185 [hep-ex]. 

42. ATLAS Collaboration, Updated luminosity 
determination in pp collisions at ^/s = 7 TeV using the 
ATLAS detector, ATLAS-CONF-201 1-011 (2011) , 
http : //cdsweb . cern . ch/record/1334563/. 

43. P. Nason, A new method for combining NLO QCD with 
shower Monte Carlo algorithms, JHEP 11 (2004) 040, 
arXiv:hep-ph/0409146. 

44. S. Frixione, P. Nason, and C. Oleari, Matching NLO 
QCD computations with Parton Shower simulations: 
the POWHEC method, JHEP 11 (2007) 070, 
arXiv: 0709. 2092 [hep-ph]. 

45. B. P. Kersevan and E. Richter-Was, The Monte Carlo 
Event Generator AcerMC version 3.5 with interfaces to 
PYTHIA 6.4, HERWIG 6.5 and ARIADNE 4.I, 
hep-ph/0405247. 

46. ATLAS Collaboration, Expected Performance of the 
ATLAS Experiment - Detector, Trigger and Physics, 
pp. 870-883, (2009) , arXiv:0901.0512 [hep-ex]. 

47. P. Z. Skands, Tuning Monte Carlo Generators: The 
Perugia Tunes, Phys.Rev. D82 (2010) 074018, 
arXiv : 1005 . 3457 [hep-ph] . 

48. M. Botje et al., The PDF4LHC Working Group Interim 
Recommendations, arXiv: 1101 .0538 [hep-ph]. 

49. Tevatron Electroweak Working Group Collaboration, 
Combination of CDF and DO results on the mass of 
the top quark using up to 5.8 fb-1 of data, (2011) , 
arXiv: 1107.5255 [hep-ex]. 

50. F. A. Berends, H. Kuijf, B. Tausk, and W. T. Giele, On 
the production of a W and jets at hadron colliders. 
Nuclear Physics B357 (1991) 32 - 64. 

51. T. Gleisberg et al.. Event generation with SHERPA 
1.1, JHEP 02 (2009) 007, arXiv: 0811. 4622 [hep-ph]. 

52. L. Lyons, D. Gibaut, and P. Clifford, How to combine 
correlated estimates of a single physical quantity, Nucl. 
Instrum. Meth. A270 (1988) 110. 

53. A. Valassi, Combining correlated measurements of 
several different physical quantities, Nucl. Instrum. 
Meth. A500 (2003) 391-405. 

54. J. H. Kiihn and G. Rodrigo, Charge asymmetries of top 
quarks at hadron colliders revisited, arXiv: 1109.6830 
[hep-ph] . 

55. J. Aguilar-Saavedra and M. Perez- Victoria, Simple 
models for the top asymmetry: Constraints and 
predictions, JHEP 1109 (2011) 097, arXiv: 1107.0841 
[hep-ph] . 

56. J. Aguilar-Saavedra, Single top quark production at 
LHC with anomalous Wtb couplings, Nucl. Phys. B804 
(2008) 160-192, arXiv: 0803. 3810 [hep-ph]. 

57. ATLAS Collaboration, Search for same-sign top-quark 
production and fourth- generation down-type quarks in 



16 



pp collisions at sqrt(s) = 7 TeV with the ATLAS 
detector, arXiv: 1202.5520 [hep-ex]. 

58. W. Hollik and D. Pagani, The electroweak contribution 
to the top quark forward- backward asymmetry at the 
Tevatron, arXiv: 1107.2606 [hep-ph] . 

59. V. Ahrens, A. Ferroglia, M. Neubert, B. D. Pecjak, and 
L. L. Yang, The top-pair forward-backward asymmetry 
beyond NLO, Phys. Rev. D84 (2011) 074004, 

arXiv: 1106.6051 [hep-ph]. 



17 

The ATLAS Collaboration 

G. Aad''^ B. Abbott"o, J. Abdallah", A.A. Abdelalim^s , A. Abdesselani"^ O. Abdinoyio, B. Abi^", 

M. Abolins^^, O.S. AbouZcidi^^^ H. Abramowiczi^a^ H. Abreu^", E. Accrbi*^®^'^*^, B.S. Acharyai63a,i63b^ 

L. Adamczyk^'^, D.L. Adams^"*, T.N. Addy^^ J. Adelman^'^*, M. Aderholz^^, S. Adomcit^^, P. Adragna^^, 

T. Adyei28^ g^ Aefsky22, J.A. Aguilar-Saavedrai^sb.a^ M. Aharrouche^o, S.P. Ahlen^i, F. Ahles^^, A. Ahmad"7^ 

M. Ahsan^'o, G. Aiellii32a,i32b^ rp^ Akdogan^**^, T.P.A. Akesson^^ G. Akimotoi^*, A.V. Akimov ^^, 

A. Akiyama^e, M.S. Alam^, M.A. Alam^^ J. Albert^^®, S. Albrand5^ M. Aleksa^^, I.N. Aleksandrov^^, 

F. Alessandria^^^, C. Alexa^^^, G. Alexander^^^ G. Alexandl•e■*^ T. Alexopoulos^, M. Alhroob^o, M. Alievl^ 

G. Alimonti**«^ J. Alison"9, M. Aliyeyio, P.P. AUport^^ S.E. Allwood-Spiers^^^ J. Almond^i, A. Aloisioioi^'i°i'^, 
R. Aloni™, A. Alonso^®, B. Alvarez Gonzalez^^ M.G. Alviggiioia,ioib^ j^ Amako6^ P. Amaral^^, C. Amelung^^, 
V.V. Ammosov^^'^, A. Amorim^^^'^''', G. Amoros^^^, N. Amram^^^, C. Anastopoulos^^, L.S. Ancu^*', 

N. Andari"*, T. Andeen^^, C.F. Anders^^, G. Anders^®^, K.J. Anderson^o, A. Andreazza^^^'^^'^, V. Andrei^^^, 

M-L. Andrieux^^, X.S. Anduaga*'^, A. Angerami^'^, F. Anghinolfi^^ , A. Anisenkov^^'', N. Anjos^^^'', A. Annovi*'', 

A. Antonaki^, M. Antonelli"'^, A. Antonov^^, J. Antos"^^^ F. Anulli"^'', S. Aoun*^^^ l^ Apcrio Bella-*, 

R. Apollc"^'^ G. Arabidzc^^, I. Araccna"2, y. Arai^^, A.T.H. Arce^^, J.P. Archambault^s, S. Arfaouii^^^ 

J-F. Arguin", E. Arik**^^'*, M. Ariki*'", A.J. Armbruster^^ , O. Arnaez^o, C. Arnault""*, A. Artamonov^*, 

G. Artonii3ia,i3ib^ O Arutinov^o, S. Asai*^*, R. Asfandiyarov*^*, S. Ask^^, B. Asmani*5a,i45b^ L. Asquith^ 

K. Assamagan^*, A. Astbury*''^, A. Astvatsatourov^^^ g Aubert*, E. Auge"*, K. Augsten*^'', 

M. Aurousseau"*^, G. Avolio*^^ R. Avramidou^ D. Axcn^^^^ c. Ay^*, G. Azuelos^^,^^ y. Azumai54^ 

M.A. Baak29, G. BaccaglioniSS^, C. Baccii^Sa.isab^ ^.M. Bach**, H. Bachacoui^s, K. Bachas^^, G. Bachy^^, 

M. Backes*^, M. Backhaus^o, E. Badescu^^^^ p. Bagnaia"ia,i3ib^ g^ Bahinipati^, Y. Bai^^^, D.C. Bailey*^'^, 

T. Baini57^ j.t. Baincsi^s, O.K. Baker*^*, M.D. Baker^*, S. Baker™, E. Banas^^, P. Banerjee^^^ Sw. Banerjee*^*, 

D. Banfi29, A. Bangertl*^ V. Bansal*^®, H.S. Bansil*^, L. Barak*™, S.P. Baranov^^^ A. Barashkou^*^ 

A. Barbaro Galtieri*", T. Barber''^ E.L. Barberio^^ D. Barberis^O'^'^'"", M. Barbero^^, D.Y. Bardin^^, 

T. Barillari^*, M. Barisonzi*™, T. Barklow*'^^ N. Barlow^'^, B.M. Barnett*^^, R.M. Barnett**, A. Baroncelli^^Sa, 

G. Barone*^, A.J. Barr**'', F. Barreiro™, J. Barreiro Guimaraes da Costa^^, P. Barrillon***, R. Bartoldus**^^ 

A.E. Barton™, V. Bartsch*"*, R.L. Bates^^^ l^ Batkova*43a^ j ji Batley^^, A. Battaglia*^ M. Battistin^^, 

F. Bauer"5, H.S. Bawa**2'^ S. Beale^^ B. Beare*^^ T. Beau^'^, P.H. Beauchemin*60^ ^ Beccherle^°^, 

P. Bechtle^o, H.P. Beck*^, S. Becker^^, M. Beckingham"^^ K.H. Becks*™, A.J. Beddall*^^, A. Beddall*^^ 

S. Bedikian*^*, V.A. Bednyakov^*, C.P. Bee^^ M. Begel^^, S. Behar Harpaz*^*, P.K. Behera^^, M. Beimforde^s, 

C. Belanger-Champagne^*, P.J. Bell*^, W.H. Bellas, G. Bella*^^^ L. Bellagamba*^^, F. Bellina^^, M. Bellomo^^, 
A. Belloni^'^, O. Beloborodova*"^./, K. Bclotskiy9^ O. Beltramello^^, S. Ben Ami*'"^*, O. Benary*^^ 

D. Benchekroun*^***, C. Benchouk^^^ y[^ Bcndel^°, N. Benekos*^*, Y. Bcnhanimou*'''^, E. Benhar Noccioli*^, 
J.A. Benitez Garcia*^^*^, D.P. Benjamin**, M. Bcnoit***, J.R. Bensinger^^, K. Benslama*^^, S. Bentvelsen*°*, 
D. Berge^^, E. Bergeaas Kuutmann**, N. Berger"*, F. Berghaus*^*, E. Berglund*°*, J. Beringer**, P. Bernat™, 
R. Bernhard*^, C. Bernius^'*, T. Berry™, C. Bertella®^ A. Bertin*^^'*^^, F. Bcrtinelli^a, F. Bcrtolucci*2i^'i2ib^ 
M.I. Besana®^'*'*^'", N. Besson*3^ S. Bcthke9^ W. Bhimji*^ R.M. Bianchi^^, M. Bianco^*^'^**^, O. Biebel^^ 
S.P. Bieniek™, K. Bierwagcn^*, J. Bicsiada**, M. Biglictti*33a, H. Bilokon*^ M. Bindi*^'^'*^^, S. Binct***, 

A. Bingul*^^ C. Bini*3ia'i3ib^ c. Biscarat*™, U. Bitenc*^, K.M. Black^*, R.E. Blair^, J.-B. Blanchard*35, 

G. Blanchot29, T. Blazek**^^^ c. Blockcr22, J. Blocki3^ A. Blondel*^ W. Blum^o, U. Blumenschein'^*, 

G.J. Bobbink*04, V.B. Bobrovnikov*^^ S.S. Bocchetta™, A. Bocci*^ C.R. Boddy**'', M. Boehler**, J. Boek*™, 
N. Boelacrt^s, S. Boser™, J.A. Bogaerts^^, A. Bogdanchikov*^^, A. Bogouch^^'*, C. Bohm**^^, V. Boisvert™, 
T. Bold^^ V. Boldea^sa, N.M. Bolnet*^^, M. Bona'^"^ y.G. Bondarenko^^ M. Bondioli*^^^ M. Boonekamp*35, 
G. Boorman™, C.N. Booth*^^, S. Bordoni'^^ C. Borer*^, A. Borisov*^'', G. Borissov™, I. Borjanovic*^^, 
M. Borri**, S. Borroni^^, K. BosiO"*, D. Boschcrini*^'^, M. Bosnian**, H. Boterenbrood*°'', D. Botterill*^*, 
J. Bouchami^2^ j Boudreau*^^, E.V. Bouhova-Thacker™, D. Boumediene'^'^, C. Bourdarios**^, N. Bousson®^, 
A. Boveia^", J. Boyd^^, I.R. Boyko^'*, N.I. Bozhko*^^, I. Bozovic-Jelisavcic*^^, J. Bracinik*'^, A. Braem^^, 
P. Branchini*33a, G.W. Brandenburg^'^, A. Brandt^ G. Brandt**'^, O. Brandt^'*, U. Bratzler*^^, B. Brau^^ 
J.E. Brau**3, H.M. Braun*™, B. Brelier*^^ J. Bremer^^, R. Brenner*^^, S. Bressler*™, D. Breton**", 
D. Britton^^^ P.M. Brochu^^ I. Brock^o, R. Brock^^ T.J. Brodbeck™, E. Brodet*^^, F. Broggi^*'^, 
C. Bromberg^^, J. Bronner^^, G. Brooijmans'^*, W.K. Brooks"^*'^, G. Brown^*, H. Brown^, 
P.A. Bruckman de Renstrom'^*, D. Bruncko**'^'^, R. Bruneliere*®, S. Brunet^°, A. Bruni*^^, G. Bruni*^^, 



18 

M. Bruschii^^, T. Buanes", Q. Buat^^ F. Bucci''^ J. Buchanan"^, N.J. Buchanan^, P. Buchholz^o, 

R.M. Buckingham"^, A.G. Buckley''^, S.I. Buda^^^, LA. Budagov^", B. Budick^o^, V. Biischer^o, L. Bugge"^ 

0. Bulckov^^ M. Bunsc-*2, T. Buran"^, H. Burckhart^^, S. Burdin'^^^ T. Burgcss^^, S. Burkc^^s^ g^ Busato^^, 
P. BusseyS^^ C.P. Buszcllol'^^ F. Butin^^, B. Butlcr"2, J.M. Butler^i, CM. Buttar^^^ J.M. Butterworth™ , 
W. Buttinger^^, S. Cabrera Urbanise, D. Caforio^^'^'i^'', O. Cakir^^, P. Calafiura", G. Caldcrini^^ 

P. Calfayan^^, R. Calkinsl°^ L.P. Caloba^^^, R. Caloii^i^asib^ j^ Calvet^^, S. Calvet^^, R. Camacho Toro^^ 
P. Camarri"2a,i32b^ ^^ Cainbiaghi"^^'"^'^, D. Cameron"^^ l.M. Caminada", S. Campana^^, M. Campanelli^^^ 
V. Canalei°i^'i°i'", F. Canelli^o^s, A. Canepa^ssa^ j^ Cantero™, L. Capassoi°i^'i°i'", M.D.M. Capeans Garrido^^, 

1. Caprini^sa, ^^ Caprini^s^, D. Capriotti^^, M. Capua^e^'^eb^ ^ Caputo^°, C. Caramarcu^^, R. Cardarelli"2a^ 
T. Garlics, G. Carlinoioi^, L. Garminati**^^'^**'^, B. Garon^^, S. Garon^o^^ G.D. Garrillo Montoyai''\ 

A.A. Garter^4^ J.R. Garter^^, J. Garvalhoi^a^,'!^ d. Gasadei^^^ M.P. Gasado", M. Gascellai^i^^izib^ 

G. Gaso5°^'^°^'*, A.M. Gastaneda Hernandezi^\ E. Gastaneda-Miranda^^i, V. Gastillo Gimenezi^e^ 

N.F. Gastroi23a^ G. Gataldi^^'', F. Gataneo^^, A. Gatinaccio^^, J.R. Gatmore^^, A. Gattai^^, G. Gattami32a,i32b^ 

S. Gaughron*^, D. Gauzi*^3a,i63c^ p^ Gavalleri^^, D. Gavalli^^^, M. Gavalli-Sforza", V. Gavasinmi2ia.i2ib^ 

F. Gcradini"3a,i33b^ ^ g^ Gerqucira^sb, A. Gerri^^, L. Gerrito^^^ F. Gcrutti^^, S.A. Getini^'^, F. Cevenmi^^^'''^^^^ , 

A. Ghafaqi34a^ d. Ghakrabortyl°^ K. Ghan^, B. Ghapleau^^^ j p^ Ghapman^^, J.W. Ghapman^e, E. Ghareyre^^, 

D.G. Gharltoni^ V. Ghavda^^, G.A. Ghavez Barajas^^, S. GheathamS^, S. Ghekanov^, S.V. Ghckulacyi'^Sa^ 

G.A. Ghelkov'^4^ M.A. Ghclstowskaio^, G. Ghen^^ H. Ghcn^^, S. Ghen32<=, T. Ghcn32c, ^ chcni^\ S. Ghcng^^^^ 

A. Gheplakov'^4^ y p^ Ghepurnov^^, R. Gherkaoui El Mourslii34c^ y. Ghernyatin^*, E. Ghcu^, S.L. Ghcungi", 
L. Ghcvalicr"^ G. Chiciari^"^^'^"^^ , L. Ghikovani^i% J.T. Ghilders^^, A. Ghilingarov™, G. Ghiodini7l^ 

M.V. Ghizhov64, Q Choudalakis^o, S. Ghouridoui^e, LA. Ghristidi^^^ a. Ghristov^s, D. Ghromck-Burckhart^s, 
M.L. Ghu^^o, J. Ghudobai24, G. Giapctti"i^^"i'^, K. Giba^^, A.K. Giftci^^, R. Giftci3% D. Ginca^^, V. Gindro^^^ 
M.D. Giobotaru^s^ G. Giocca^^^, A. Giocio", M. GiriUi^^ M. Gittcrio^®'', M. Giubancan^s^, A. Glark'*^ 
P.J. Glark^s, W. Glelandi22, J.G. Glemens82, B. Glement^s, G. Glement"5a,i45b^ j^^^ Gliffli^s, Y. Goadou*^^^ 
M. Gobali63a,i63c^ j^ Goccaroi^i, J. Gochran^s, P. Goe"^ J.G. Gogan"2^ j. Goggeshalli64^ E. Gognerasl^^ 
J. Golas^ A.P. Golijnio^^ N.J. GoUinsi^ G. Gollins-Tooth^s, J. GoUot^s, G. Golon^^^ P. Gonde Muifioi^sa^ 

E. Goniavitis"^, M.G. Gonidi", M. Gonsonniio^^ V. Gonsorti4^ S. Gonstantinescu^Sa, G. Gonta"^^'"^'^, 

F. Gonventiioi'^'% J. Gook^^, M. Gooke^^, B.D. Gooper™, A.M. Gooper-Sarkar"^, K. Gopic^^, T. Gornclisscn^^^ 
M. Gorradi^^'', F. Gorriveau^^J', A. Gortes-Gonzalez464^ G. Gortiana9^ G. Gosta*^**'', M.J. Gostal^^ 

D. Gostanzo^^s, T. Gostin^o, D. Gote^^, R. Goura Torres^^^, L. Gourneyea^^^s, G. Gowan^^ G. Gowden^^, 
B.E. Gox*\ K. Granmerio^ F. Gresciolii2ia,i2ib^ ^^ Gristinziani^o, G. Grosetti^^a.aeb^ ^ Q,^pj7ia,7ib^ 

S. Grepe-Renaudin'^^, G.-M. Guciuc^''^, G. Guenca Almcnar^'^'', T. Guhadar Donszelmann"^ , M. Guratolo'''^, 
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A.D. Hershenhorn^^^, G. Herten*®, R. Hertenberger^^, L. Hervas^^, N.P. Hessey^^'', E. Higon-Rodriguez*^^, 
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J.L. Holzbauer^^, Y. Homma'^^ T.M. Hong"^ L. Hooft van Huysduynen^o^ T. Horazdovsky^^e^ C. Horn"^ 

S. Horner^s, J-Y. Hostachy^^^ S. Houi^o, M.A. Houlden^^ A. Hoummada"4a^ j. Howarth^i, D.F. Howell"^, 
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M.S. Kim2, P.C Kim"^ S.H. Kinii^s, N. Kimurai^s, 0. Kindts, B.T. King™, M. King^e, R.S.B. King"^ 
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Y. Kulchitsky89, S. Kuleshov^ib^ C. Kummer^^, M. Kuna^^ N. Kundu"^ J. Kunkle"^ A. Kupco^^^^ 
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C. Le Maneri57^ e. Le Menedeu^, C. LebeF^ T. LeCompte^ F. Ledroit-Guillon55, H. Lee^o^, J.S.H. Lcg"^^ 
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